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Abstract 



We calculate the static potential, the drag force and the jet quenching parameter 
in strongly coupled anisotropic M = 4 super Yang-Mills plasma. We find that the jet 
quenching is in general enhanced in presence of anisotropy compared to the isotropic 
case and that its value depends strongly on the direction of the moving quark and 
the direction along which the momentum broadening occurs. The jet quenching is 
strongly enhanced for a quark moving along the anisotropic direction and momentum 
broadening happens along the transverse one. The parameter gets lower for a quark 
moving along the transverse direction and the momentum broadening considered 
along the anisotropic one. Finally, a weaker enhancement is observed when the 
quark moves in the transverse plane and the broadening occurs on the same plane. 
The drag force for quark motion parallel to the anisotropy is always enhanced. For 
motion in the transverse space the drag force is enhanced compared to the isotropic 
case only for quarks having velocity above a critical value. Below this critical value 
the force is decreased. Moreover, the drag force along the anisotropic direction is 
always stronger than the force in the transverse space. The diffusion time follows 
exactly the inverse relations of the drag forces. 

The static potential is decreased and stronger decrease observed for quark-antiquark 
pair aligned along the anisotropic direction than the transverse one. We finally 
comment on our results and elaborate on their similarities and differences with the 
weakly coupled plasmas. 



1 Introduction 



Last years there is a lot of effort to understand the heavy-ion collisions and the 
quark-gluon plasma (QGP). The Relativistic Heavy Ion Collider (RHIC) findings 
and the following analysis suggests that the QGP is a strongly coupled fluid [1] and 
therefore the perturbative methods of Quantum Chromodynamics (QCD) are in 
general not appropriate for describing it. On the other hand there is some progress 
in Lattice field theory, (see for example [2] ) but further progress is very difficult since 
in QGP we need to study real-time phenomena. Moreover, a promising approach to 
study these phenomena is developed with the use of methods within gauge/gravity 
duality [3] , where an up to date review is in reference [4] . Using the AdS / CFT it is 
possible to study several observables and properties of the dual QGP working in the 
strongly coupled regime. At the moment the studies are mostly in the qualitative 
level, but still the information that can be extracted is significant. 

Although the exact gravity dual of the QCD is not known, and the theories 
and their dual backgrounds used for QGP calculations usually have different char- 
acteristics and properties than QCD (although some of them disappear in high 
temperatures), there are some important results that appear to have some kind of 
universality among the different theories. Relativistic hydrodynamics describe well 
the QGP [5] and one of the most known results obtained so far is for the prediction 
of the ratio of shear viscosity over entropy density [6]. Moreover several methods 
have been developed for calculations of the jet quenching, the drag force and the 
relaxation time, the static potential and the quarkonia etc. in the dual QGP, which 
accompanied with some arguments for normalizing the results so to lead to more 
sensible comparisons with the results obtained from the experiments. 

By now there is a large number of papers which study these quantities in the 
phase where the plasma has become already isotropic and at equilibrium. However 
the plasma after its creation and for a short period of time is anisotropic both in 
momentum and coordinate space. For the RHIC energies the elliptic flow of the mat- 
ter created is described quite well by models that assume that the hydrodynamical 
behavior is applicable at r < lfm. The conformal viscous hydrodynamics predicts 
times r ~ 2fm but the estimation depends strongly on the initial conditions (eg. 
Color Glass Condensate (CGC) or Glauber) and details of plasma hadronization. On 
the other hand, by considering the collision of two sheets of energy in strongly cou- 
pled M = 4 sYM in [7] has been predicted a thermalization time of order ~ 0.3 fm. 
Therefore, the current estimations of thermalization time vary significantly. 

In this paper we initiate the study of several observables in a dual anisotropic 
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strongly coupled QGP. Such anisotropy that we study here, can be referred to the 
momentum space and is caused due to locally anisotropic hydrodynamic expansion 
of the plasma. We point out that the question we answer accurately in this paper 
is how the observables are modified in the dual M = 4 sYM plasma in presence 
of anisotropy. Whether or not our results apply to the observed anisotropic QGP, 
depends mainly on how well the initial isotropic theory, of which our theory here 
is a consistent deformation, describes the QGP. There are several indications that 
at least qualitatively the predictions in the isotropic case capture characteristics of 
the real QGP quite well, so our analysis here as well could capture properties of 
the anisotropic plasma, and we find that our results are indeed consistent to what 
is expected. Apart from that, it is very interesting on its own to see how several 
quantities in context of AdS/CFT are modified in presence of such an anisotropy 
we consider here. 

Recently there is important progress in the anisotropic gauge/gravity dualities. 
In [8] has been found an anisotropic dual geometry with a naked singularity. This 
geometry was used in [9] to study electromagnetic signatures of the dual plasma. 
In [10] the anisotropic supergravity solutions that were found are regular on and 
outside the horizon. The geometry was used in [11] to find that the longitudinal 
shear viscosity over entropy bound is violated. The strong coupling isotropization 
of a large number of anisotropic initial states in the absence of external sources was 
also studied recently in [12]. In this paper we continue and extend the studies on 
the anisotropic gauge/gravity dualities. 

Regarding the QGP anisotropics, it is known that the existence of them is im- 
portant for the QGP evolution. For example, a spatial anisotropy which occurs due 
to the fact that in the heavy ion collision the nuclei have finite area and usually col- 
lide off-center causes the hydrodynamic elliptic flow. The flow is associated to the 
elliptic flow parameter 112 which is defined as the anisotropy of particle production 
with respect to the reaction plane and is a way to measure how the system responds 
to the initial spatial anisotropy. The elliptic flow is generated because the pressure 
gradient along the impact vector direction on the transverse plane is larger than the 
one in the transverse direction and the nuclear medium expands preferably along 
the impact vector direction. Moreover the interaction between the particles, leads 
to a momentum anisotropy distribution on the reaction plane reflecting the excited 
medium to the above spatial geometry. The measurement of the elliptic flow pro- 
vides information for the thermalization times and can be used constrain the ratio 
of shear viscosity over entropy density. 

Here we are mainly interested for the anisotropy which created by the rapid ex- 
pansion of the plasma along the longitudinal beam axis at the earliest times after 
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the collision. The longitudinal pressure to the beam axis is lower than the trans- 
verse one and the momenta of the partons along the beam direction are lower than 
the ones in the transverse space. By considering boosted hadrons with the same 
velocity the effective temperature increase with the mass of the hadron species [13]. 
These momentum distribution anisotropies cause plasma color instabilities which 
are responsible for the isotropization short time and process of the QGP [14]. 

All the different anisotropies mentioned can occur at the same times. A way 
to isolate the anisotropy we are interested on is to think the colliding nuclei as 
having infinite transverse area, or that the collisions are completely central. After 
the collisions the partons are produced at the formation time where the partonic 
momentum distribution can be supposed to be isotropic. Then a rapid longitudinal 
expansion of the plasma along the beam line occurs. During this process the lon- 
gitudinal expansion rate is larger than the parton interaction rate, and the plasma 
along the longitudinal direction is much colder than the one in the transverse di- 
rection. At this stage the pressure along the longitudinal and transverse directions 
satisfy Pl < Pt and the corresponding momenta (p 2 L ) < (p^) in the local rest frame. 
At the time r = r iso the interaction rate becomes equal to the expansion rate the 
plasma reaches the isotropic phase where the hydrodynamic analysis can be done. 
This momentum anisotropic plasma have chromo-Weibel instability, which believed 
to play important role on the isotropization process. 

In this paper motivated by the experimental as well as theoretical studies, we 
initiate the study several observables in a dual anisotropic plasma. We use a high 
temperature limit of a static, regular IIB supergravity solution dual to a spatially 
anisotropic finite temperature M = 4 super Yang-Mills (sYM) plasma [10]. The ge- 
ometry characterized by an anisotropic parameter where its limit to zero is smooth 
and gives the isotropic undeformed finite temperature M = 4 sYM. This is ex- 
pected since the supergravity solution can be seen as a deformation of the original 
solution, where the anisotropy can be though as introduced either by a non-zero 
number density of dissolved branes that do not extend to boundary and therefore 
do not add new degrees of freedom, or resulting from a #-term which depends on 
the anisotropic direction. In the dual gravity the particular #-term corresponds to 
an axion depending on the anisotropic direction. 

Using this background we start by calculating the static potential between a pair 
of heavy probe quarks QQ. We find how the critical length of the pair distance and 
the values of the static potential depend on the anisotropic parameter along the 
different directions in the plasma. We compare these results to the isotropic case 
and then compare with models that study the static potential in the weak coupling 
regime. Moreover, another motivation for these calculations is that this study might 
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be useful to extract qualitative results for the quarkonia in the anisotropic plasma. 

Continuing we calculate the drag force and the diffusion time of a heavy quark 
moving along different directions of the anisotropic plasma. We derive the analytical 
results for quarks moving along the transverse and longitudinal directions, compare 
them each other and to the isotropic case. 

Then we study the jet quenching parameter. It's bounds can be measured in 
the QGP by the radiative energy loss and the parameter itself can be though as a 
property of the strongly coupled medium. In our 4-dim anisotropic plasma we have 
three different choices for the transverse momentum broadening. The energetic par- 
ton moves along one of the transverse directions and the momentum broadening 
happens along the anisotropic direction. The second is when the parton moves par- 
allel to the anisotropic direction and the momentum broadening considered in the 
transverse direction. For the last one the parton moves along the transverse to the 
anisotropy directions and the momentum broadening is calculated along the other 
transverse direction. In this case although only the transverse directions are con- 
sidered, the dependence of the radial metric element on the anisotropic parameter, 
modifies sightly the result compared to the one in the undeformed theory. In gen- 
eral we find enhancement of the jet quenching in presence of anisotropy. Our results 
again compared to the results obtained from other models. Finally, we also discuss 
the difficulties that arise in our model when we try to give a more precise quantita- 
tive prediction using different comparison schemes. These difficulties appear because 
we are working on a small anisotropy over temperature limit. 

In this paper we have tried to present some of the analytic calculations clearly in 
the Appendices in order to improve the readability of the main text. The structure 
of the paper is as follows. In the second section we present the background and the 
theory we use and how its parameters are related to the parameters of models using 
anisotropic momentum distribution functions. In the following section we investigate 
the static potential in the anisotropic dual plasma. This section is supported by the 
Appendix A, where the generic gravity dual orthogonal Wilson loop calculations are 
presented. In section 4, the drag force and the quark relaxation time is studied. 
The appendix B supports this section where the analytical calculations for the drag 
force are presented for any general background. In section 5, we study the jet 
quenching. Here we present in the main text the jet quenching calculation and 
the approximations done to derive the result for a complete generic background, 
since the calculation is very interesting and certain approximations done in the 
derivation need to be tested in our background. Then we derive the results in 
the anisotropic background and comment on them. In the next section we make 
an attempt to provide more quantitative predictions of our results using different 
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comparison schemes. We finalize with the discussion session where we also collect 
our results. 

2 The model 

2.1 The dual geometry 

The anisotropic background we use is a deformed version of the M = 4 finite tem- 
perature sYM [10]. The deformation parameter in the field theory is introduced 
by a ^-parameter term depending on the anisotropic direction. It turns out that 
9 = 2^7107X3, where x% is the anisotropic gauge theory coordinate and ud 7 can be 
thought as the density of D7-branes homogeneously distributed along the anisotropic 
direction. The 9 angle is related to the axion of the type IIB supergravity through 
the complexified coupling constant of the M = 4 sYM. Therefore in the gravity dual 
background the anisotropic deformation can be seen as inserted due to existence of 
axion term depending on the anisotropic direction, or as due to the backreaction 
of the D7-branes. These branes do not add new degrees of freedom to the theory 
since they do not touch the boundary. They are wrapped on the internal space and 
the transverse directions to the anisotropy, therefore creating the anisotropy on the 
deformed AdS geometry. 

In the string frame the background is given by 

ds 2 = i (^-FB dx 2 + dx\ + dx\ + Udx\ + ^-^j + Z dQ 2 SB . (1) 

X = ax 3 , (f) = (j)(u) , (2) 

where a is the anisotropic parameter with units of inverse length, <p is the dilaton, 
X is the axion depending linearly on £3 coordinate. The anisotropic direction is 
considered to be the X3 and the functions J 7 , £>, H depend on the radial coordinate 
u and the parameter a. The background has also a RR five form but it is not 
important for our purposes. The analytical form of the functions can be found for 
small anisotropy compared to the temperature or sufficiently high temperatures, 
T 3> a. It is enough to consider the expansions of the fields up to second order in a 
around the black D3-brane solution: 

F(u) = l-^j + a 2 F 2 {u) + 0{a A ) (3) 

B{u) = 1 + a 2 B 2 {u) + 0{a A ) , (4) 
%(u) = e-* w , where <j>{u) = a 2 2 (n) + C(a 4 ) . (5) 

Note that only even powers can appear because of the symmetry z —z. By 
applying asymptotic AdS boundary conditions and requiring Ti to vanish at the 
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horizon u = Uh, the Einstein equation can be solved giving: 
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We can find the temperature evaluating the following expression at the horizon 
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and solving for the horizon position Uh we get 
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As expected the isotropic limit a — >■ reproduce the well know results of the isotropic 
black D3-brane solution. The parameters of our background now are the tempera- 
ture T and the anisotropy a. The energy density per unit volume can be calculated 
from 
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where gL4/i is the area element of the hypersurface £ = const, u = Uh- In our case 
the result reads 

+ 0(a 4 ). (10) 

The energy and pressures can be found from the expectation value of the stress 
tensor, where the element (T 00 ) is the energy, (T n ) = (T 22 ) = P Xl x 2 ='■ P± denote the 
pressure along the x± or x 2 direction or the transverse plane and (T 33 ) = P X3 =: Py 
is the pressure along the anisotropic direction, calling it also longitudinal one. The 
analytic expressions read 
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Therefore for high temperatures the pressure of the plasma along the anisotropic 
direction is always lower than the one in the other two directions 



P < P 

1 X3 ^ 1 X\X2 



(12) 
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Figure 1: Anisotropy in the momentum space for £ > 0. 



2.2 Relation of background parameters to anisotropic mo- 
mentum distribution function 

In weakly coupled anisotropic plasmas a usual technique used for the study of the 
observables is to consider an anisotropic phase space distribution function. The 
degrees of freedom in weakly coupled plasmas are split to soft modes that carry 
momenta of order gyhiT and hard ones that carry moment of order T. The hard 
modes are particles that have an anisotropic phase space distribution function. 

To make connection of this function with our parameters we can consider a 
kinematic example where the accelerated beams with nucleons collide along the X3 
anisotropic direction. This is the beam-axis direction and along this direction the 
system expands initially. The plasma created can be seen as having a distribution 
function f(t, x, p) which can be taken homogeneous in position space but anisotropic 
in momentum space. The anisotropic distribution function that can be written 
as [15] 



where the vector n = (0, 0, 1) is the unit vector along the anisotropic direction 
and the parameter £ plays the role of the anisotropic parameter. This distribution 
represents a stretched or contracted version of the isotropic case since one direction 
in the momentum space is rescaled. For —1 < £ < the distribution is stretched 
along the anisotropic direction while for £ > the distribution is contracted in the 
anisotropic direction as in Figure 1. 

The parameter £ is related to the average particle momenta in the transverse 
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p T = p — n(p ■ n) and longitudinal pl = p • n to anisotropy directions by 



where the factor of 2 appears in the denominator due to the number of the trans- 
verse directions. Therefore we see that the anisotropic distribution can be obtained 
from the isotropic one by removing or adding particles that have large momentum 
component along the anisotropic direction. The anisotropic plasma created after 
the heavy ion collisions correspond to values £ > 0, and this is the region that we 
work here. 



In order to find a relation between the anisotropic parameter £ and the parameter 
a of our supergravity background, we have to express the first one in terms of the 
pressures of the system. This is done by using the distribution (13) and calculating 
the pressures through the stress energy tensor components. A new parameter A 
can be introduced which measures the degree of momentum anisotropy through 
pressures and defined as 

A:=^-1 = %L-1. (15) 
After some algebra A can be related to £ [16] as 

A=i({-3) + {((l + {)5^-l)", (16) 
which in the small A and consequently small £ limit, the expression becomes 

limA = t£ + 0(£ 2 ) (17) 



and in the large £ limit 



lim A = h + 0(^) . 



For our background in the high temperature limit using (11) we get for A 
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In the range of T > a =^ A < 1 we can relate the two anisotropic parameters as 

5a 2 „ , 
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where £ has to be positive and very small. The equation (20) provides the basic con- 
nection between the parameters of our supergravity background and the anisotropic 
momentum distribution functions (13) considered in several field theory models. 



Moreover if f iso represents an ideal gas momentum distribution and £ is small 
enough, the anisotropic parameter is related to shear viscosity over entropy density 
and to the proper time of the plasma. For one dimensional Bjorken expansion the 
analytic relation is [17] 

1 TS 

where the anisotropy increases with the expansion rate. In the following sections we 
calculate several physical observables, we explain our results and try to qualitatively 
compare our results with experimental data and the weakly coupled plasma models. 



3 QQ Static Potential in the anisotropic M = 4 
plasma 

In this section we study the static potential in the finite temperature anisotropic 
dual plasma. In order to do so we are using the analytic equations (91) and (92) 
derived in the Appendix A. We will compare the static potential to the isotropic 
case as well as the potentials along the different anisotropic directions. 

We use the usual ansatz for the string world-sheet choosing the static gauge and 
the a dependence on the radial direction: 

xq = r and x p = a, u = u(a) , (22) 
where x p = X\ =: x± or x p = x 3 =: x\\ . (23) 

In the first case we align the QQ pair along the direction of x\, which is equivalent 
aligning it along x-i- Then the pair is placed along the anisotropic x^ direction. 
Finally in order to compare with the isotropic finite temperature M = 4 sYM theory 
we set a = and calculate the static potential, where the analysis is equivalent to 
the configuration considered in [18]. These are the three different static potentials 
we study and compare each other. 

Initially we fix the temperature and the anisotropy parameter and consider the 
static potential of quark pairs with different separation lengths. By applying the 
formulas derived in the Appendix A we find that the potential for the pairs aligned 
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Figure 2: The static potential close to L c for 
QQ pairs aligned along the anisotropic and 
the transverse direction with respect to L. 
Their relation is | Vn \ < \V±\. The correspond- 
ing critical lengths are L c | < L c ±. Settings: 
blue dotdashed line-Vy , red solid line-V^ and 
T = 3, a = 0.3T. 
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Figure 3: The ratios of the static poten- 
tials for pairs aligned along different direc- 
tion. All the fractions reduce as the pair 
distance increase and approach the criti- 
cal length. Settings: V\\ /Vj_-black dotdashed 
color, Vj|/Vi so -blue dashed line, Vj_/Vi so -solid 
red line and T = 3, a = 0.35T. 
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Figure 4: Vji for two values of the anisotropic 
parameter. The static potential decreases for 
increasing the anisotropy. Settings: a\ = 
0.5T, a 2 = 0.01T and T = 3. 



Figure 5: The VTi and V± where the string 
world-sheet turning point is fixed to uq /uh = 
0.5 and a increases. The o^0 point is where 
the lines cross, and as a increases the lines di- 
verge. By comparing the two final points of 
the curves which correspond to a ~ 0.5T we 
see that as a increases the background geom- 
etry and the potentials along the parallel di- 
rection is affected more than transverse one. 
Settings: as in Figure 2 and T = 3. 



10 



along the anisotropic direction tends to be slightly weaker compared to transverse 
alignment as can be seen in Figures 2 and 3 l : 

V\\ < V± < V iso for — 7^ and constant. (24) 

We also find that the difference between the potentials in two directions gets bigger 
as the distance of the quarks approach the critical length L c . Moreover, the critical 
length is reduced in presence of anisotropy as 

L c || < L c j_ < L c iso for -^0 and constant. (25) 

In order to study the static potential dependence on the anisotropy parameter we 
keep constant the temperature and vary the anisotropic parameter. Increase of 
anisotropy leads to decrease (Figure 4) of the absolute value of the static potential 
in any direction compared to the isotropic case and therefore to decrease of the 
critical length 

a/ ^V l± \ =>L cl± \ . (26) 

A more detailed analysis shows that as the anisotropy increases, the deviations of the 
anisotropic static potentials along the different directions increase (Figure 5). The 
anisotropy affects the configuration along the anisotropic direction stronger than 
the transverse one. This can be seen in Figure 5, where we fix the temperature, the 
ratio uo/uh and so consequently the length L of the Wilson loop with respect to the 
horizon position, and we increase the parameter a. Notice that now along the rr-axis, 
same lengths L along different directions, correspond to different values of anisotropy 
but the first and the last points in the two curves correspond to the same value of 
the anisotropy parameter. We see that V\\ = V± for a — > 0, and by increasing a the 
potentials and the pair distances along the different directions deviate. This happens 
because the different strength on the modifications of the background geometry along 
the different directions is reflected to the static potential. 

Finally, by keeping the anisotropy constant and increase the temperature the 
static potential in any direction gets lower absolute value and the critical length 
reduces (Figure 6). This is not unexpected since our dual theory is a smooth defor- 
mation of M = 4 sYM where the same behavior has been observed [18]. 

1 We should comment here that in some of the plots we choose the anisotropic parameter a of 
order 0.3T or so. Although smaller values give the same behavior, their effects are not visible in 
the plots clearly. 
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Figure 6: The form of the potential in anisotropic background is the same as the unde- 
formed N = 4 sYM one. This is not unexpected since the anisotropy introduced by a 
continuous deformation to W = 4 sYM. Settings: T = 1 and a = 0.001. 

3.1 Comments and comparison with other models 

To interpret our findings one could think that lower pressure in the parallel direction, 
which results having less energetic gluons leads to the further screening of the po- 
tential. However this seems not to be the case. The enhancement of the anisotropy 
while keeping constant temperature, results to increase of P Xl x 2 an d simultaneous 
decrease of P X:j as can be seen from (11), or equivalently leads to contraction of 
the momentum in the anisotropic direction, by removing for example energetic glu- 
ons in that direction. At the same time momentum extension happens along the 
transverse direction and increase of the relevant pressure. However, here we find 
that in both directions the absolute value of the static potential is decreased. An 
other observation that indicates that the static potential is not related directly to 
the pressures is that at least to a 2 the pressures are modified by the same amounts 
from the initial one. If the force modification between the quarks would depend 
mainly on the pressure then the magnitude of modification should be almost equal 
in parallel and transverse directions. This does not happen and it is not expected 
to happen judging from the differences in the metric in the two directions. The 
decrease of the static potential is more likely to be related to the increase of energy 
or entropy density of the system as found in equations (10) and (11). 

Notice that at some point in our analysis we have compared the potential in 
anisotropic background to the isotropic one. When comparing the potentials along 
the parallel and transverse directions in the anisotropic background for fixed a/T the 
horizon of the black hole remains at the same position and only the corresponding 
geometry that the string worldsheets extend is modified. While comparing to the 
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isotropic case the position of the black hole horizon differs but the comparison can be 
made by identifying the temperatures in the isotropic and anisotropic backgrounds. 

Our theory does not have dynamical degrees of freedom in the fundamental 
representation. Inclusion of quarks in isotropic M = 4 sYM with additional D7 
flavor branes beyond the probe limit is expected to lead to screening of the potential 
at least at scales close to the string breaking one [19,20]. When the anisotropy 
introduced in a background that contains already dynamical degrees of freedom is 
not clear from our analysis whether or not the potential will be further reduced. This 
will depend on the flavors dependence of the anisotropy. If for example the density 
of flavors turn out to depend on the anisotropy and it is reduced for increasing 
anisotropy then the static potential could increase since screening will be weaker 
and the decrease due to gluons found here will be an antagonizing behavior. On 
the other hand if the density of the flavors will be independent of the anisotropy, 
then the potential even in presence of flavors will be reduced as we have found 
here. Additionally, the relation between the static potential along the different 
directions might get modified after the inclusion of flavors, since their dependence 
on the anisotropic parameter and therefore their contribution to screening can not 
be a priori predicted. Therefore in this paper we have found the behavior of the 
static potential in presence of anisotropy and matter in adjoint representation, but 
a further analysis is needed in presence of unquenched dynamical quarks. 

It is interesting to report in this section some findings of the static potential in 
the weak coupling regime. In [21] it has been found stronger attraction for quark 
pairs aligned along the direction of the anisotropy than for transverse alignment and 
that the potential generally gets enhanced in presence of anisotropy. These results 
are valid on distance scales on the order of the inverse of the Debye mass (while 
the relevant plots are for L ~ (0.5, 2)m~j D 1 ) and for small anisotropic parameter. In 
this regime a naive extrapolation shows that already in our configuration we have 
entered in the deconfined phase. Moreover the authors find that as the separation 
length of the pair gets reduced the difference between the potential in the transverse 
and parallel directions is reduced. The static potential results in our strong coupling 
theory without flavors differ with the ones observed for the weak coupling plasmas. 
However there exist several differences between the models apart of course for the 
main difference that these findings are in weak coupling. 

Before we move on to the drag force analysis we summarize only some of the 
findings. In this section we have found that the anisotropy in our model results to 
the decrease of the static potential and affects more significantly the potential for 
a pair aligned along the anisotropic direction. The static potentials along different 
directions follow the relation (24) and increase of the anisotropy leads to decrease 
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of the static potential in any direction. The critical lengths £ c ||,± depend on the 
anisotropy and for increasing anisotropy they reduce (25). 

4 Drag force on the heavy quarks in the anisotropic 
plasma 

The drag force in an anisotropic plasma depends on the direction of the motion 
of the heavy external probe quark. We consider a heavy quark moving along the 
anisotropic direction and then along the transverse direction with a fixed velocity 
v, in an infinite volume of gluon plasma and fixed finite temperature T [22]. As 
usual the force we measure here is the one that needs to be imposed on the quark in 
order to keep it moving with the constant velocity v. This force can be generated for 
example from a constant electric field, which will balance the backward drag force 
on the quark generated from the interaction of quark with gluons and quarks in the 
plasma. 

The velocity v is bounded in order the drag force calculation to be valid. To 
make sure that the quark indeed loses energy and the generated force from the 
plasma is a backward drag force, it should move well above the subthermal velocities 
v 'subthermai ~ \/T/Mq, where Mq is the mass of the heavy quark. So the low bound 
in the velocity is w > v su bthermai- Additionally there is an upper bound which 
should be imposed in order for our calculation to be valid. This is required because 
for example to keep very large constant speed v through the plasma a very large 
electric field is needed which after a critical value would produce qq pairs. There are 
different ways to calculate the maximum value of the velocity: either by consider 
the DBI action of a D7-brane that represents the heavy quark and taking a reality 
condition or by comparing the deceleration forces from drag and vacuum radiation 
and imposing the first one to be the dominant. The maximum value of the velocity 
in our theory remains approximately similar to the undeformed M = 4 sYM and is: 



In order to describe the trailing string for a motion along the x p := x\\ t ± directions 
we use the radial gauge choice: 




(27) 



where VXT/Mq < I. 



X = T, 



u = a. 



x p = vr + f(u) , 



(28) 
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where in the other directions the world-sheet is localized. The drag force derivation 
can be done generically for any background. We derive the force in the Appendix 
B. The final equations we obtain and we use in this section are the equations (97) 
and (99). It is interesting that the drag force depend on the time metric element, on 
the metric element along the direction where the quark moves and only indirectly 
to the radial metric element through (97). 

The results obtained in this section are analytical and we have also check them 
numerically. The radial position uq, where the numerator and denominator change 
sign simultaneously in (96) in order to keep f(u)' real, can be written as: 

"o|| = «oi|| + a 2 uo2\\ , u Q ± = u oi ± + a 2 u 02 ± , (29) 

where w i an d u 02 can in principle depend on all the parameters of the problem, but 
not on the anisotropic parameter a. They are considered for parallel and transverse 
motion with respect to the anisotropy parameter. Plugging them in the equation 
(97) for parallel and transverse motion correspondingly and solving each of them 
separately and for different orders, we obtain: 

(l-v 2 ) 1/4 

W01|| = "Oil = -j, = U iso , (30) 

VT^ 2 (1 + VT^ 2 ) + (7v 2 - 5) log (1 + VT^?) 
U ° 2 " " 48.3T3 (1 - ^ ' (31) 



VT^P (l + VT^ 2 ) + (4w 2 - 5) log (l + VT^) 

U ° 2± ~ 48^(1-^ • (32) 

Notice that the Woi is equal to the isotropic result as expected and the only difference 
in uq2 expressions for different directions is a numerical factor. The drag forces using 
equation (99) are given by: 

Fdrag, \\ nT 2 V 



drag,. 

~7T 





(33) 



(34) 



v 2 



2(1 - v 2 ) - (5 - Av 2 ) log (1 + Vl^) 



!) (1 + ^1^2) (1 _ ,,2)3/2 



The 0th order term in anisotropy is equal to the drag force in the isotropic unde- 
formed background as it should be. The anisotropic correction does not depend on 
the temperature in any direction and it is consistent with the dimensional analysis. 
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Figure 7: The drag force dependence on the 
anisotropic parameter a. The velocity is cho- 
sen close to 1, v ~ 0.98. The exact value of 
the velocity as far as is in this region, does not 
affect qualitatively the results. All the ratios 
in this plot are greater than one but for v < 
0.909 the Fdrag,±/ Fdrag,iso is lower than unit 
and reduces as the anisotropy increases. Set- 
tings: black dotdashed ]m&F drag \\/F dragt ±, 
blue dashed Hne-F drag ^/F drag:iso , red solid 
\m.e-F drag ^/F dragMo and T = 1. 



Figure 8: The drag force dependence on the 
quark velocity v for fixed a and T. The only 
quantity lower that 1 is the F drag ±/ F dragtiso 
for v < 0.909. Settings: as in Figure 7 and 
a = 0.1 and T = 1. 




Figure 9: The drag force ratios dependence 
on the temperature T. We choose v = 0.98, 
where all ratios are above unit. As the tem- 
perature increases the ratios decrease. Set- 
tings: as in Figure 7 and a = 0.1. 



Figure 10: The drag force ratios dependence 
on the temperature T. We choose v = 0.8 
where the ratio F drag>± / F drag:iso < 1. As the 
temperature increases the deformed forces 
approach to the undeformed one. Settings: 
as in Figure 7 and a = 0.1. 



The a 2 term in the F drag \\ is always negative for any value of v indicating en- 
hanced drag force for motion along the anisotropic direction. The corresponding 
term in F drag ^± is positive for velocities smaller than v < v c ~ 0.909 indicating de- 



1(3 



crease of the drag force in this region, compared to the undeformed force. Above this 
critical value the drag force even in the transverse direction is enhanced, compared 
to the isotropic one as can be seen from (34). It is also worth noticing that the value 
of the critical velocity is independent of the temperature and anisotropy. 

Moreover the drag force in the anisotropic direction is always greater than the 
one in the transverse direction and their ratio reads: 



The correction term is always positive and depends on the temperature as expected 
from the dimensional analysis. 

For better understanding of our findings we plot the following results. In the 
Figure 7 are presented the drag forces along parallel and transverse directions 
normalized with the undeformed isotropic result Fd ra g,iso, as well as the fraction 
Fdrag, \\ / F^ag, ± i depending on the anisotropy. We see that increase of the anisotropy 
leads to increase of the deviation of all the drag forces. The strongest dependence 
on anisotropy is observed for quarks moving along the anisotropic direction as was 
also observed for the static potential. 

In the Figure 8 we plot the same quantities depending on the velocity of the quark 
for constant a and T. We observe that increase of the velocity of the probe quark 
leads to increase of the drag force. Where again for high enough velocities the drag 
force for motion along the anisotropic direction is affected stronger by the anisotropy. 
From the plot it is clear that for v < 0.909 the drag force in an isotropic background 
is stronger than the force on the quark moving along the transverse direction to 
anisotropy. For higher velocities all the drag forces are enhanced compared to the 
isotropic case. 

In Figures 9, 10 we plot the drag force ratios for fixed anisotropy and velocities. 
We observe the differences for v < v c and v > v c . Common in both plots is that as 
the temperature increases with respect to anisotropy the drag forces approach the 
isotropic ones as expected by construction of our background. 

To summarize qualitatively our results: 

Fdrag, \\ Fdrag,iso and Fd r ag,\\ ^* Fdrag, _L > 

Fdrag,± > F dr ag,iso for v > v c , while below this velocity 

Fdrag,!. ^ Fdrag,iso- 




drag, | 



1 + a 



2 (2 - v 2 ) log (1 + VT^) 
8tt 2 T 2 (1 - v 2 ) 



(35) 
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4.1 The drag coefficient and the diffusion time 



The drag coefficient is denned as 



d P -44, M Q V fO(l\ 

Tt =-nnP, with P=^=j- (36) 



Therefore the diffusion time tjj is given by: 



Td,\\,X = = ~ n 2 ^ 



n D,\\,± Fdrag,\\,± 



V 



and corresponds to the time where the initial momentum is reduced by e _1 factor. 
The equation (37) indicates that the time inequalities between the different direc- 
tions will be just inverted compared to the drag force results. For example the ratio 
of the diffusion time for motion along parallel and transverse direction is given by: 



^,11 I , (2-i> a )bg(l + vT=^) 

^~ l ~ a 8vr 2 T 2 (1 - v 2 ) ' (38) 

meaning that the diffusion time along the anisotropic direction is lower compared 
to the one in the isotropic direction. The other two ratios are 2 : 

r D ,\\ a 2 v / r^(l + v / T^) + (l + ^ 2 )log(l + v / r^) _ 

— 1 — — . — > 



T DMo T 2 24tt 2 (1 - V 



2^ 



Tn iso T 2 247T 2 (1-^) " 1 ] 



The diffusion time td,± is longer than the isotropic one for v < v c , while for v > v c 
it is shorter as expected from the drag force results. In Figures 11, 12, 13 we plot 
the ratios of diffusion times with respect to velocity of the quarks for three different 
values of a. Increase of the velocity of the moving quark or increase of the anisotropy 
makes the ratios TD,\\/TD,(i S o,±) to increasingly diverge from the unit. 



To briefly summarize, we find that for very high velocities the diffusion time 
in presence of anisotropy is reduced compared to the isotropic theory. Hence the 
diffusion process is faster in presence of the anisotropy for high velocities. We also 
find that the anisotropy affects more the diffusion time in the anisotropic direction 

2 Notice that here we are considering as negligible the modifications to the mass of the quark 
from the thermal medium. As we discuss in section 6, the medium induced corrections to the mass 
arc similar in any direction, but are expected to be slightly modified compared to the isotropic 
case. 
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Figure 11: The ratio of the diffusion time 
in the transverse direction of the anisotropic 
holographic model to the diffusion time of the 
isotropic theory. For v < 0.909 the ratio is 
above 1. Settings: red solid line-a = 0.1, blue 
dashed line-a = 0.2, black dotdashed line- 
a = 0.3 and T = 1. 
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Figure 12: The ratio of diffusion time for 
motion along the anisotropic direction to the 
diffusion time of the isotropic theory for dif- 
ferent values of the anisotropic parameter. 
Settings as in Figure 11. 




Figure 13: The ratio of diffusion times along 
the anisotropic direction to the transverse one 
is always lower than the unit. For higher ve- 
locities the ratio diverges increasingly from 
the unit. Settings as in Figure 11. 



that in the transverse one. More particularly the relation between the different 
diffusion times are: 

t\\ < r iso and r\\ < r± , (41) 
t± < r iso for v > v c , while below this velocity r± > r iso . (42) 
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5 Jet Quenching in the anisotropic strongly cou- 
pled Af = 4 sYM plasma 



In this section we calculate the anisotropy effects on the jet quenching parameter. Its 
bounds can be measured in the QGP by the radiative energy loss and the parameter 
itself can be though as a property of the strongly coupled medium. The jet quenching 
is generated when the momentum of an energetic parton changes while it moves in 
the medium. The interaction with the medium result the parton to radiate gluons 
which is the reason of the transverse momentum broadening. The jet quenching 
is defined as the fraction of the mean transverse momentum obtained by the hard 
parton in the medium over the distance it has traveled. 

In the field theory the transverse momentum broadening can be calculated using 
a Wilson line in the adjoint representation along a light cone direction. In the gravity 
dual description the jet quenching can be calculated from the minimal surface of a 
world-sheet which ends on an orthogonal Wilson loop lying along two light-like lines. 
These two long parallel lines of the Wilson loop, with length say L_, are related 
to the partons moving at relativistic velocities and are taken to be much more 
larger that the other two sides of the loop with length related to the transverse 
momentum of the radiated gluons. 

The Wilson loop we calculate in supergravity side is in fundamental represen- 
tation but in the planar limit the expectation value of the adjoint Wilson loop is 
related to the Wilson loop in the fundamental representation as 

(W A (C)) = {W F {C)f , (43) 

because trAdj = t r F un d- The LHS of the above equation is related to the jet quenching 
parameter as [23] 

(W A (C)) « exp"^ L ' L - . (44) 

To calculate the corresponding Wilson loop we go to the light-cone coordinates by 
the coordinate transformation \f 7 lx ± = t ±x p , where x p is chosen to be x\\ or x±, as 
in (85) . The generic metric (83) becomes 

ds 2 = G (dx 2 + + dx 2 _) + G + -dx + dx- + Ga^ p )dx 2 + G uu du 2 (45) 

G — = - (Goo + G pp ) , G + - = Goo — G pp 

Taking advantage of the condition L_ ^> L^, where Lk is the length of the string 
in the k direction, we assume that the string worldsheet is translational invariant 
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along x k . This simplifies things significantly. We present here the full calculation 
of the jet quenching because under certain approximations we can arrive in an 
analytic expression valid for any background. We numerically check our analytical 
expressions with the exact expressions and we find that our approximations are 
correctly considered. 

The ansatz for the string configuration is 

£_ = r, Xk — cr, u = u(a) (46) 
x+, x^p are constant , (47) 

which represents a Wilson loop extending along the Xk direction and lying at a 
constant x + ,x^ p . The indices p, k here denote a chosen direction. The action then 
reads 

Lu Li- 

S=—=; day/G—{G w u* + G kk ) =: — / day/D (48) 
27ra' Jo 27ra' J 

and the analysis has some similarities to the QQ action analysis presented Appendix 
A. Using Hamiltonian formalism we obtain 

H = . (49) 



The Hamiltonian is a constant of motion, and we set it equal to c. We can solve for 
the v! and get 

u n = (G kk G - c*)G kk (5Q) 

where at the turning point the string satisfies 

G kk = 0, or G kk G^ = c\ or G^ = . (51) 

Usually the interesting equation in this case is the last one which is satisfied for 
u = uq = Uh- The short length of the string is then given by 



t = / * du \j (G,J G '-c> )Gik ■ (52) 

Since we are interested in the small L k length and the integral goes from the bound- 
ary to the turning point, the constant c must be very small. Therefore, we expand 
our formula, and at the end we will check the validity of the expansion numerically 
to justify it. The constant of motion turn out to be given by the analytic expression 




1 G, 



-i 

-3 \ 
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Before we substitute to the action we need to eliminate the infinity in the action 
that appears due to the bounds of the integral. We use the mass subtraction scheme 
where we subtract the two straight string world sheets described by: x_ = r, u = a. 
The self energy reads: 

-- / duy/G—G uu . (54) 

JO 



So 



2ira 

The total action subtracted the divergences is equal to: 



run 

- / dusj G UU G — 
* Jo 



S - So = — / duy/G w G— ( K \ - — fc \ 2 - 1 1 . (55) 



27ra' ./n " \ V G — - c 

After some algebra the normalized action takes the simple form 



Aira 

where expressed in metric elements 



1 



S - S = - Lk c , (56) 




(57) 



We have checked numerically that in our case the small c approximation is valid 
where the equations (55) and (57) take the same values for small values of c. The 
jet quenching parameter for an energetic parton moving along the p direction while 
the broadening happens along the k direction is given by: 




Uk) = — 7 / Tc:\\n— ' ^ 



where to normalize correctly we have taken into account the equation (43). The 
equation (58) is very useful, since can be applied directly to any background that 
satisfies the approximations we have made here. 



For our 4-dim anisotropic plasma we have three different choices for the trans- 
verse momentum broadening. The first one, q\\(±), is for the energetic parton mov- 
ing parallel to the anisotropic direction and the momentum broadening occurring 
along the transverse direction. The second one, q±n\), is for an energetic parton 
moving along X\ or x<i direction and the momentum broadening happens along the 
anisotropic £3 direction. The last one, q±(±), is for a parton moving along the 
transverse to the anisotropy directions and the momentum broadening is considered 
along the other parallel direction. In this case although only the directions along the 
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transverse space are involved, the dependence of the radial direction metric element 
on the anisotropic parameter, modify although sightly the result compared to the 
isotropic theory. Therefore the configurations we consider are: 



q 


Xp 




Energetic parton moves along 


Momentum broadening along 


9J-CII) 


x± 


x \\ 


X± 


x \\ 


3||(±) 


X\\ 


x ± 


X\\ 


x± 


9-J-C-L) 


X±,i 


X±,2 


X±,l 


X±,2 



Where the x^± are defined with respect to the anisotropy direction. 



Let us start with the case where the quark moves along the anisotropic direction 
and the momentum broadening occurs along the transverse one. Then 

G pp = G 33 and G kk = G n . (59) 

We find that the jet quenching in presence of anisotropy is enhanced with respect to 
the isotropic case. Stronger anisotropy leads to further enhancement of the fraction 
q\\(±)/qo, where go : = <jiso- Notice that normalization with g_i_(_L) instead of go, leads 
to almost identical results. In Figure 14 we show some representative results. The 
form of the function is linear and particularly: 

1 + 0.122^ (60) 

in a very good approximation, where the numerical factor does not depended strongly 
on the other parameters of our background. Comparing to the other jet quench- 
ing parameters g_L(||)> Q±(±) we wm see that the anisotropy affects stronger the jet 
quenching for the quark moving along the anisotropic direction than the other ones. 

To examine the q±(\\) we need to set 

Gpp = G u and G kk = G 33 (61) 

and substitute to (58). Some representative results are shown in the Figure 15. 
Again we observe enhancement of the jet quenching, but this time is much more 
weaker than the g||(_i_). The form of the function is no more linear, but increase of 
anisotropy still leads to increase of the jet quenching. 

Finally, we look at q±t±) where we need to choose 

Gpp = G u and G kk = G22 ■ (62) 

We show that a weak enhancement of the jet quenching observed as can be also seen 
in Figure 16. Although the quark motion and the momentum broadening happen 
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Figure 14: The jet quenching for a heavy 
quark moving along the anisotropic direction 
9|i(_L) as a function of a/T. The results are 
normalized with the isotropic jet quenching. 
Notice that normalization with gives 
almost identical results. Settings: T = 5. 
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Figure 16: Quark motion and broadening in 
transverse space:^^). Settings: T = 5. 
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Figure 15: The momentum broadening q±(\\) 
along the anisotropic direction for a heavy 
quark moving in the transverse space. Set- 
tings: T = 5. 
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Figure 17: Comparing the jet quenching pa- 
rameters. The <7||(j_) is much more bigger than 
the other two as expected. The q±(\\), <7_i_(_i_) 
are closer and are plotted here. Settings: 
Blue solid line-g^(||), red dashed line-g_|_(j_) , 
T = 5. 



along the transverse directions the jet quenching gets modified due to the fact that 
the time and the radial metric elements are dependent on the anisotropy parameter. 
We compare qj_(±) with g_i_(||) i n Figure 17 where we observe important differences in 
their magnitudes. 

There are some common features that we observe for all the jet quenching pa- 
rameters. We find that when anisotropy is present the jet quenching is enhanced. 
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Increasing the anisotropy parameter leads to increase of the enhancement. The 
presence of anisotropy affects mostly the momentum broadening of a quark moving 
along the anisotropic direction. This is something we have seen in the static poten- 
tial and the drag force. Alignment of the QQ pair along the anisotropic direction, 
as well as motion of the quark with the trailing string along this direction led to 
maximum modifications on the relevant results due to anisotropies. 

To summarize our findings: We find that the jet quenching is in generally en- 
hanced in presence of anisotropy compared to the isotropic case and that its value 
depends strongly on the direction of the moving quark and the direction which the 
momentum broadening occurs. More particularly 

q\\(±) > q±(\\) > q±(x) > Qiso , (63) 

which in our conventions translated as: the jet quenching is stronger enhanced for a 
quark moving along the anisotropic direction and momentum broadening happens 
along the transverse one. The parameter gets lower for a quark moving along the 
transverse direction and the momentum broadening considered along the anisotropic 
one. A very weak enhancement is observed when the quark moves in the transverse 
plane and the broadening happens on the same plane. 



5.1 Comparison with other results 

It has been shown in many different studies that the jet quenching in an anisotropic 
plasma is increasing in presence of anisotropy. 

It has been found using kinetic theory that in the leading-log order the jet quench- 
ing of a heavy quark moving along one of the transverse directions to anisotropy is 
larger when the momentum broadening is along the anisotropic axis, which in this 
case coincide with the beam axis, compared to when the broadening occurs in the 
other transverse direction [24]. In our conventions these results can be written as 
> Q±(±)- This is in agreement with our results as can be seen in Figure 17. 

In [25] the jet quenching was calculated in an unstable non-Abelian weakly cou- 
pled SU(2) plasma. Using numerical simulations taking into account hard elastic 
collisions and soft interactions mediated by classical Yang-Mills fields, and a separa- 
tion scale between them, the authors found that the fields develop unstable modes 
which lead to q±(\\) > Q±(±) in our notation. 

In [26] the jet quenching is estimated in leading logarithmic approximation by the 
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broadening of the massless quark interacting via gluon exchange. The energetic hard 
quark considered propagates in one of the transverse directions to the anisotropy 
and the momentum broadening is estimated in the transverse to the motion plane, 
which include the anisotropic direction. The jet quenching in the anisotropic plasma 
was found to be enhanced with respect to the isotropic one. 

It is important to notice here that our jet quenching results are in agreement 
with STAR findings, eg. [27]. This might be an additional indication that the QGP 
in LHC and RHIC is not in equilibrium. 



6 Attempts for more Quantitative Comparison 

In this section we make an attempt to predict the effect of the anisotropy on the ob- 
servables in a more quantitative way. To do that we need to use sensible comparison 
schemes and give realistic values to the parameters of the model. There are several 
differences between our anisotropic deformed theory and the QCD but a compari- 
son between strongly coupled plasmas in the two theories could be done under some 
logical normalization schemes. For example one motivation for these comparisons 
is that there are indications from lattice QCD calculations which show that lattice 
QCD thermodynamics are in good approximations with conformal theories for some 
finite range of temperatures T > 2T C . 

For the heavy test quarks we used we can choose the charm and bottom quarks. 
Their effective masses are difficult to be determined precisely in the thermal medium, 
but the most representative values are: M c = 1.5GeV and = 4.8GeV. We 
can include the medium induced correction to these masses, by specifying their 
dependence on the temperature. As we mention in the Appendix A, the mass 
of the static quark is represented by the static straight string stretched along the 
radial radiation starting from the boundary of the space and reaching the black 
hole horizon. However, this string is infinite and that is the reason that is used 
to normalize the infinite static potential. Therefore in the UV, a regulator should 
be introduced to make the result finite. The value of the regulator has to be fixed 
by matching the Mq(T = 0) to the physical quark mass. Notice that the medium 
induced corrections to the masses are not affected of the direction of the anisotropy. 
On the other hand, they are very slightly modified compared to the isotropic case, 
since the horizon position is modified. 

To compare the observables in the theories we need to fix appropriately the theory 
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parameters. Our model is not a confining and the static potential does not include 
the linear term. However, by comparing the renormalized charge or the static force 
of the conformal theory with the lattice data for relatively small separation lengths 
of the quark pair, it is possible to find that a well estimated value in order to agree 
for the t' Hooft coupling is A = 5.5. The next parameter we need to determine for 
comparison is the choice of temperature in our anisotropic deformed theory and the 
QCD. The degrees of freedom in two theories as well as the field content are very 
different. For example in our case we do not have flavors and we have very large 
number of the color branes. To qualitatively compare the two theories we normalize 
the quantities firstly according to a fixed energy density scheme. By approximating 
the QCD plasma clS db free gas the energy density is 

€ qgp * ^ NcNf) T* QGP ~ 11.2 T* GP , where N c = N f = 3 , (64) 

and the factor iV c 2 - 1 + N c N f is the degrees of freedom of SU(N C ) QCD with N f 
flavors. By counting the degrees of freedom in M = 4 sYM can be found that there 
are approximately 2.7 times more than the ones in QCD above the confinement 
phase transition, in our case ~ 45. Taking advantage of the fact that our theory 
is deformed with a very small deformation parameter a and equating the energy 
densities in these two theories we obtain the relation [28] 

TsYM,anisot = 2.7^'Tqcd ■ (65) 
The next parameter we need to fix is the anisotropy. From (20) we find that 

* = *5!f ( 66 ) 
5 

and therefore we are in the region of small £. If the anisotropic parameter £ is small 
then it is also related to the shear viscosity of the plasma, where in one-dimensional 
boost invariant expansion governed by Navier-Stokes evolution [16,30,31], the result 
given from (21) 

5 = (67) 

1 TS 

where the rj/s is the viscosity to entropy ratio, and r _1 defines the expansion rate, 
since r is the proper time. The dimensionless parameter Tr determines the mag- 
nitude of the anisotropics and in the strong coupling in the RHIC and LHC initial 
conditions respectively this can be approximated to 

Tr ~ 0.35 for RHIC conditions (68) 
Tt ~ 0.43 for LHC conditions , (69) 
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where we have assumed tq ~ 0.2/m, T = 350MeV and To ~ 0.1/m, T = 850MeV^ 
respectively, since at LHC energies the initial time is expected to be small. To get a 
better picture of the above assumptions and the qualitative picture we fix the total 
entropy and observe the behavior of Tr. By assuming that the entropy scales as 
s oc T 3 and V oc r , with S = s V we can estimate the corresponding times at 
different conditions. For example, by considering To = 250 MeV to correspond to 
tq — 1, then for temperature To = 350 MeV the corresponding time is tq ~ 0.36. 
Therefore, increase of temperature leads to decrease of To which is in agreement 
with the assumptions made to get (69). By choosing in the strong coupling limit 
that rj/s ~ 0.1 as a representative value, we get for the initial conditions £ ~ 2.3 
and £ ~ 2.8, where the smaller momentum space anisotropies correspond to the 
LHC initial conditions. Notice that reduction of the number of the flavor degrees of 
freedom, result also lower values for To since the temperature increases. 

For RHIC energies the typical formation times cover a range of 0.3/m to less 
than Ifm for temperatures T = 250 — 350 MeV, estimation made by fixing the total 
entropy to reproduce the measured rapidity density of hadrons at a given centrality, 
e.g., at impact parameter b ~ Ifm [29]. Additionally, the hydro dynamical models 
normally require thermalization times Ttherm in the range of 0.6 — Ifm in order to 
reproduce the magnitude of elliptic flow which is observed at RHIC. We need to 
choose a representative value for £ for T = 250 MeV and by considering r ~ 0.6/m, 
then Tr ~ 0.76 and 

U ^ 1.30 . (70) 

This is the value we choose as a representative one for RHIC energy. Notice that 
the choice of the value depends strongly on the values we chose for the viscosity 
over entropy ratio and for the the proper time. Since the exact values for both 
of these quantities are quite uncertain our numerical value should be taken with 
caution. For the viscosity over entropy ratio the value that is predicted in the 
AdS/CFT is around the well known r]/s ~ 1/4tt [6]. For RHIC conditions it has 
been estimated for a hot gluon plasma for SU(3) pure gauge model using lattice 
QCD simulation a range of values 0.1-0.4 [32]. More recently for the SU(3) gauge 
theory an estimate of rj/s > 0.134(33) is predicted at T = 1.165T C [2]. Therefore, a 
choice of t]/s ~ 0.1 as representative is sensible. The proper time value has chosen 
to have this value in order to be larger than formation time and close to the lower 
bound of the thermalization time so to estimate the anisotropy effects before the 
thermalization. 

The thermalization time at the LHC is estimated lower than the RHIC, at least to 
0.1/m [29]. Then by choosing a representative value r ~ 0.5 fm and T = 450MeV, 
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relevant for LHC conditions we get Tr = 1.14 and 



Cl ^ 0.87 . (71) 

Therefore in order to compare between our anisotropic theory and the QCD, the 
last thing we need to specify is the relation between their proper times. One first 
approximation would be to consider proper times approximately equal. An other 
option is to fix the total entropy of motion since in our model the anisotropy pa- 
rameter is very small. This assumption seems to be logical and in agreement with 
the discussion for heavy quarks in [29]. In that case the formation times between 
two theories with different degrees of freedom follow the relation 

Tqcd \ 3 ( d-o.f-QCD 



T(W=4 sYM — I 7f I 7 T 0QCD ■ (72) 

V J A^=4 sYM J \ a -0-J-M=4: sYM J 

Then we get for the formation times 

tqm=a sYM — 2.7" 1/4 t qc*d (73) 

and for the Tr 

(Tr )^ =4 sYM = 2.7-^ 2 (Tr ) Q cD ■ (74) 

We see that both the formation time and the Tr quantity is lower in M = 4 
(anisotropic) sYM. This implies for £ close to the initial conditions 

^N=i sYM ^ V2?7 £ Q cd , (75) 

while the numerical values are for the LHC and RHIC initial conditions we chose 
are translated to £r ^=4 sYM — 4.60 and £l m=4 sYM — 3.78. 



We can assume that the proper time relation (72) between the two theories 
carry on beyond the formation time, for any proper early time we would like to 
compare: ri^ =4 s ym = ^-^^^Tiqcd for our comparison purposes. For example this 
can be though as measuring the time with respect to the formation time in the 
two theories. In that case £r a SYM — 2.14 for the representative RHIC value (70) 
and £l a sYM — 1-43 for the representative LHC value (71). Therefore from (66) 
the resulting values for the supergravity anisotropy parameter are (a/T)^ ~ 5.81 
and (a/T)x ~ 4.75. For these values our current background can not be used for 
predictions since we are in T > a approximation. 

Notice if we would fix the entropy density instead of the energy density the 
relation between the formation times by making the above assumption of the con- 
stant total entropy would obviously be r ^ =4: s ym — toqcd- This would lead to 
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(Tr )Af=4 sYM = 2.7~ 1/3 (Tt ) QC d and consequently to £^=4 sY m - 2.7 1/3 £q C d- In 
this case to the values of £ parameters are very close to the ones found above, so 
qualitatively for our purposes there is no change. 

Using the direct scheme where the temperatures of the two theories are identified, 
the isotropization time would be equal in the two theories since in TV = 4 sYM it 
was found that Ti SO ~ 0.7/T [33]. Then the values of £ near the isotropization would 
coincide in the two theories and take the values (70) and (71) leading to 4.53 and 
and 3.7 respectively. The a/T values turn out to be greater than one and close to the 
values predicted above. So for our purposes the interesting outcome here is that by 
fixing different quantities in the two theories the proper times and more importantly 
the Tt quatities do not differ significantly, and all of them lead to relatively high 
a/T values if we trust the equation (21). 

To make the comparison using our model we should have had at most say a/T ~ 
0.3 which gives £, a sYM — 0.005 which correspond to extremely high temperature 
according to equation (21). Therefore, we need the full analytical background and 
not only the expansion for small a/T to make a valid comparison. Notice that 
the values that we have used here are approximate, however any sensible values for 
the proper time (and how it is modified in the comparison scheme used) and the 
viscosity over entropy ratio that could be chosen would lead to the similar outcome 
for the range of values of large a/T and anisotropy £. 

The results we have obtained in our paper show how the particular observables 
are modified in presence of anisotropy in the strong coupling limit, and the pattern 
that they follow is clear. However for completeness we consider a/T ~ 0.3 and plug 
its value to our results. Then the jet quenching modified as 

~ 1.036 , ~ 1.0009 , ~ 1.0002 . (76) 

% <?o <?o 

While the drag force, for velocity v = 0.95 is 

^11 ^1 F ± 

-f ~ 1.0035 , — ^ ~ 1.0036 , -=±- ~ 1.0001 . (77) 

" A. "iso "iso 



7 Conclusions and List of Results 



In this paper we have studied the jet quenching, the drag force and the static 
potential, in the strong coupling dual anisotropic QGP. We have seen how certain 
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quantities of the finite temperature isotropic M = 4 sYM theory are modified in 
presence of anisotropy. In the QGP plasma, similar anisotropies can be created due 
to the expansion along the longitudinal direction. We have found several interesting 
results where some of them are listed briefly below: 

• Jet quenching: In presence of anisotropy the jet quenching is enhanced com- 
pared to the isotropic case and its value depends strongly on the direction of the 
moving quark and the direction which the momentum broadening occurs. The jet 
quenching is strongly enhanced for a quark moving along the anisotropic direction 
and momentum broadening happens along the transverse plane. It reduces for a 
quark moving along the transverse direction and the momentum broadening consid- 
ered along the anisotropic one. Finally, it reduces further when the quark moves in 
the transverse plane and the broadening happens on the same plane. We can write 
our findings in a more compact form as 

> ?±(||) > <?±(J_) > Qiso ■ (78) 

Additionally the dependence of <j\\(±) on the a/T is linear in a very good approxima- 
tion with 

1 + 0.122^. (79) 

• Drag Force: In presence of anisotropy the longitudinal drag force F drag i\ is 
always enhanced. The transverse force Fd ra g,±_ is enhanced above a critical velocity 
value v ~ 0.909, while below this value is reduced. In both cases the anisotropic 
terms contributing to the anisotropy do not depend on the temperature (33), (34). 
Comparing the forces between longitudinal and transverse directions we get 

F drag ,\\ _ 2 (2-i; 2 )Log + 
F drag , ± +a 8tt 2 T2(1-^) 

which indicates that F dragt \\ > Fd rag ,± ■ 

• Diffusion time: The diffusion time behavior is in direct analogy with all 
the findings of the drag force inverted. In the transverse direction t± < Tj SO for 
v > 0.909 while this relation is inverted for velocities lower than the critical value. 
The diffusion time along the longitudinal direction is always lower that the time 
in the isotropic medium. The relation between the times in the longitudinal and 
transverse directions are t\\ < t±. The analytic relations for t\\, rj_ normalized with 
the isotropic diffusion time are given by (39), (40) and their fraction by (38). 

• Static Potential: We find that the static potential in presence of anisotropy 
becomes weaker in absolute value. The relation of the static potential of a pair 
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(80) 



aligned longitudinal to the anisotropy compared to the static potential of a pair 
aligned transverse direction is 

V\\ < V x < V iso (81) 

when the comparison is done in terms of LT keeping the anisotropy parameter and 
temperature fixed. Moreover the critical length of the string is decreased in presence 
of anisotropy as 

L c \\ < L c ± < L c iso . (82) 

The modifications of all our results due to anisotropy are stronger along the 
anisotropic direction. This happens because we have found that the geometry is 
modified more in the anisotropic direction than in transverse one, phenomenon that 
gets even stronger as the anisotropy parameter is increased. The geometry modifi- 
cations reflect to the observable results. 

It is worth noticing that for large velocities the drag force and jet quenching 
parameter have the same qualitative behavior. This qualitative agreement has been 
observed in other cases too. In [34], the fraction of drag forces and jet quenching in 
cascading plasma or charged plasma over the M = 4 sYM have been found to follow 
similar patterns. It is worthy pointing out that exact numerical agreement for the 
fractions q\\(±)/qo and F drag> \\/F can be found for velocities v ~ 0.9946, partially 
reflecting to the fact that jet quenching calculation treats ultra-relativistic quarks. 
The velocity that the fractions coincide, happens to be almost the same with the 
corresponding velocity in cascading plasmas [34]. Moreover it has been found that 
the presence of R-charges enhances the jet quenching [35]. Therefore, the inclusion 
of anisotropy and R-charges is expected to lead to further enhancement. 

The particular enhancement in jet quenching due to the anisotropy and the 
longitudinal, transverse relations we have obtained are in agreement with several 
models for weak coupling plasmas, eg. [24-26] as well as with the STAR findings [27]. 
The static potential findings are different than the ones obtained in weak coupling 
limit, but in a different setup in [21] where the potential was calculated from the 
Fourier transform of the static gluon propagator in a hard loop approximation. 
Moreover in [36] for the weakly anisotropic plasmas it was also found that the 
quarkonium binding is stronger for non- vanishing viscosity and expansion rate, result 
that is in agreement with [21]. In strong coupling limit using our model where no 
dynamical flavor degrees of freedom, we find weaker potential due to anisotropy. 
Although the inclusion of flavors beyond probe (unquenched) approximation in an 
isotropic theory is expected to lead to further screening [19,20] it is not clear a priori 
what happens in case of presence of anisotropy. We have commented on the reasons 
in the section 3.1. 
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Moreover, using some sensible comparison schemes we have tried to make a more 
'quantitative' predictions and comparison with the experimental results. However, 
the limit we have considered T a and the background is analytically known, 
corresponds to a low values of anisotropic parameter £, that seem not to be the ones 
that the observed QGP has, if we trust the equation (21). 

There are several other studies that can be performed using the AdS/CFT for 
anisotropic QGP. The quarkonium physics, where the introduction of flavor branes 
are needed in order to insert to the background information for the dynamical degrees 
of freedom and the mesons. The electromagnetic observables are also interesting to 
be studied and is believed that they contain information about the initial stages 
of the anisotropic QGP. An effort to this directions has been done in [9] using the 
background [8]. The photon and dilepton thermal production by the plasma has 
been studied in the weak coupled anisotropic plasma in [37]. Very interesting is also 
the derivation of supergravity anisotropic solutions which depend on the time and 
describe the isotropization of the plasma. 

It would be also interesting to see if and how the generic conditions for the cance- 
lation of the UV divergences in the Wilson loops with the use of the Legendre trans- 
form for isotropic backgrounds, derived in [38] are modified in case of anisotropic 
backgrounds. In the special case of the orthogonal static potential Wilson loop, it 
seems that the result of the UV divergence with the use of the Legendre transform, 
gives the same formula as the isotropic case. However, for more complicated Wilson 
loops the analysis is more involved. 

Appendix A: QQ Strings in generic weakly coupled 
backgrounds 

In this section we present the world-sheet calculation of a string in static gauge 
in weakly coupled backgrounds and their energy which corresponds to the static 
potential. The string in the anisotropic background (1), is a special case of the 
strings we examine here. The following equations have been derived in various 
forms in several in other papers too [19,39]. By writing the metric of the space as 



ds 2 = Goodr 2 + Gudx 2 + G uu du 2 , 




we choose the static gauge for the string 



x = r 



and 



Xp (X, 
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which is extended in the radial direction, so u = u(a). The x p coordinate represents 
the direction along which the pair is aligned and can be chosen to be 



x p = Xi j2 =: Xj_ or 

x p = £3 =: X|| . (85) 
Then the induced metric g a p = GMNd a X M d/3X N for our world-sheet read: 

9oo = Goo, 9n — Gp P + G ull u' 2 . (86) 
Supposing that we are working in Lorentzian signature the Nambu-Goto action is 3 

S= 2^J dadr \/- G oo(G PP + G uu u' 2 ) =: ^7 J dadr^fD . (87) 
The Hamiltonian then is equal to 

H _ GooGpp tQQ ^ 



and is a constant of motion. Setting it equal to — c we can solve for u' and get the 
turning point equation 



(GopGpp + c 2 )Gpp 

u=± \ ' (89) 

y ^ "uu 

which is solved for 

G~l = , or G pp = , or G 00 G PP = -c 2 . (90) 

The above equations, normally the last one, specify how deep the world-sheet goes 
into the bulk, and we call this value of the turning point u . 

The length 4 of the two endpoints of the string on the brane is given by 



ru ° du 2 d / ~G uu c 2 
u ' Juo V (^00 + c 2 )G pp 



3 In the case of Euclidean signature, the formulas change with a minus sign wherever the Goo 
element is. 

4 The limits of the length integral depend on where we choose as the starting point measuring 
L. When the string in the boundary extends from —L/2 to L/2, then the corresponding solution 
of the v! is positive for (0,L/2) since the turning point corresponds to L = 0. When the string 
extends from (0, L) the u' in (0, L/2) is negative. In any case the final result in the definite integral 
is (91). 
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Moreover, the energy of the string using as renormalization method the mass sub- 
traction of the two free quarks is 



2na'E 



daC 



cL 




du\j —G UU G 



GppGoo 



du\Z—G 00 G l 



(92) 



where is the possible horizon of the metric. Notice that we already used that fact 
that the world-sheet is symmetric with respect to turning point u and hence the 
RHS of the above equation are already multiplied by two. 



Using the equations (91) and (92) the static potential can be found in terms of 
the distance of the QQ pair for any gravity dual background at least numerically. 
The derivation of the analytical expressions depend on whether the integrals can be 
done analytically and if the inversion of Uq(L) is possible analytically. 



Appendix B: Drag Force on trailing string in general weakly 
coupled backgrounds 



It is possible to calculate the drag force in a generic background, for example as 
in [40], and as we have done for the static potential. We consider the background 
with the generic metric (83) and the following ansatz for the trailing string with 
radial gauge choice and motion along the x p direction: 

x = t, u = a, x p = vr + f(u) , (93) 

where in the other directions the world-sheet is localized. The x p is chosen to be xn 
or x±, like the (85) for the static potential analysis. It denotes the direction along 
which the hevy quark moves with velocity v. Moreover, the function f(u) at the 
boundary should be equal to zero in order to resemble a constant quark motion. 
The Nambu-Goto action in our case reads 

S = J du^J- (G 00 + GppV 2 ) G uu - GooGppf' 2 =: -i^J J du\f~D[ . 

(94) 

The action produces one non-trivial equation for / and since it does not depend 
explicitly on /, it has a constant of motion the canonical momentum 

n » - • < 95 > 
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The /' function must be real and to elaborate on the reality condition we need to 
solve (95) for /' to get 



y/- (Goo + G pp v 2 ) 



-GooGpp (l + GooGpp (27ra'ni) 2 ) 



(96) 



To keep /' real, both the numerator and denominator must have the same sign at 
any value u. The numerator vanishes for u = Uo satisfying 



G U u(Goo + GppV 2 ) 







(97) 



where usually the solution comes from the expression within the brackets. The 
momentum then can be expressed in terms of Uo where the denominator simplified 
significantly 



n!, 



-GooGpp 



2na' 



(9l 



U=UQ 



Since we have chosen the physical solution which describes the momentum flowing 
along the string from the boundary to the horizon, the total drag force on the string 
for motion along the x p direction is given by 



drag,x p 



ni 



Vx- 



-GooGpp 



(27T) 



(99) 



U=UQ 



Therefore in any generic background the drag force can be calculated by solving 
(97) to specify the point u and then by plugging its value to the equation (99). 



Notice: 



The results of this paper were reported in a talk, at 31 January of 2012 in the 
workshop 'Exploring QCD frontiers: from RHIC and LHC to EIC, Stellenbosch 
2012, where at the same time the presentation of the author was available online in 
the workshop webpage. While this paper was in final stage of typing for submission 
the author has received the reference [41] which has small partial overlap with our 
drag force force results, and in particular there is an agreement for small a/T limit. 
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